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Abstract 

We study the low-momentum ghost propagator Dyson-Schwinger equation (DSE) in 
Landau gauge, assuming for the truncation a constant ghost-gluon vertex, as it is ex- 
tensively done, and a simple model for a massive gluon propagator. Then, regular DSE 
solutions (the zero-momentum ghost dressing function not diverging) appear to emerge 
and we show the ghost propagator to be described by an asymptotic expression reliable 
up to the order 0{q^). That expression, depending on the gluon mass and the zero- 
momentum Taylor-scheme effective charge, is proven to fit pretty well the low-momentum 
ghost propagator obtained through big-volume lattice simulations. 

UHU-FP/ 10-023 
LPT-Orsay/10-29 

1 Introduction 



A few years ago (see for instance |T|), a vanishing gluon propagator and a diverging 
ghost dressing function at zero-momentum in Landau gauge were extensively accepted as the 
solutions for the tower of Dyson-Schwinger equations (DSE). In contrast, alternative DSE 
solutions were also predicted to give a massive gluon propagator |2l[3]. Lattice QCD (LQCD) 
estimates for those propagators appeared to be also in contradiction with a gluon propagator 
that vanishes at zero-momentum or with a ghost dressing function that diverges [IHZ]- We 
addressed this issue in two recent papers [Slin] and tried to clarify the contradiction. After 
assuming in the vanishing momentum limit a ghost dressing function behaving as F{(f') ~ 
(g^)°^ and a gluon propagator as A(g^) ~ (g^)"'^-! (or, by fohowing a notation commonly 
used, a gluon dressing function as G{q^) = g^A(g^) ~ {<f')'^°), we proved that the ghost 
propagator DSE (GPDSE) admits two types of solutions: 



• \i ap 7^ 0, the low- momentum behaviour of both gluon and ghost propagators are 
related by the condition 20^-1- ac = implying that F'^{q^)G{q^) goes to a non- 
vanishing constant when — > 0. 

• If Qi? = 0, the low-momentum leading term of the gluon propagator is constrained not 
any longer by the leading but instead by the next-to-leading one of the ghost propagator, 
and LQCD solutions indicating that F'^{q^)G{q^) when q^ ^ Q [31[5] can be pretty 
well accomodated within this case. 

^Unite Mixte de Recherche 8627 du Centre National de la Recherche Scientifique 
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In particular, the numerical study in ref. [H] of the GPDSE using a LQCD gluon input finds 
that both cases of solutions appear depending on the value of the strong coupling constant 
at the renormalization point, which is a free parameter in this exercise. Indeed, it seems to 
be by now well established that the two classes of solutions, dubbed "decoupling" (ap = 0) 
and "scaling" [ap ^ 0) may emerge from the tower of DSE [2l|3l|T0]. Such a nomenclature, 
despite being widely accepted, can be misleading. The perturbative running for the coupling 
constant renormalized in Taylor-scheme is given by F'^{q^)G{q^) and one can thus define, 
although not univocally in the IR, a coupling with it. Nevertheless, neither a scale invariance 
nor a decoupling of the IR dynamics for the theory can be inferred from the low-momentum 
behaviour of such a coupling. In particular, as will be seen, an effective charge can be properly 
defined for phenomenological purposes such that it reaches a constant at zero-momentum in 
the decoupling case. However, although not appropriate for phenomenological purposes in 
the IR domain, the Taylor-scheme coupling is a very convenient quantity in discriminating 
the kind of solutions we deal with. 

On the other hand, it was also proved in ref. [8j that, for an appropriate coupling constant 
value at the renormalization momentum, the resulting ghost dressing function (belonging 
to the decoupling class) fits very well with lattice results. It is worth pointing too that 
lattice data can be also very well accomadated within DS coupled equations in the PT- 
BFM scheme [2l[3] and within the Gribov-ZwanzigeiH approach leading in both cases to 
decoupling solutions for gluon and ghost propagators. 

Furthermore, in ref. [9], the low-momentum first correction to the constant leading be- 
haviour of the ghost dressing function for the decoupling solution was proven to be propor- 
tional to logg^ when the zero-momentum gluon propagator is constant {ac = 1), as lattice 
data seems to points to (very recentely, the authors of [H], in a different context, have also 
found a ghost propagator dressing function whith the same low- momentum behaviour). The 
proportionality factor in front of it was also proven to be written in terms of the coupling at 
the renormalization momentum and the zero-momentum values of the gluon propagator and 
ghost dressing function. The aim of this note is to go further in the low-momentum analysis 
for the ghost propagator behaviour in that decoupling case with ac = 1. With this purpose, 
a simple model for a massive gluon propagator (where the gluon mass is taken not to run 
with the momentum and to be approximated by its zero-momentum value) is applied in order 
to compute the C'(g^)-correction for the low-momentum ghost dressing function. Then, we 
prove that this low-momentum behaviour is controlled by that gluon mass and by the zero- 
momentum value of the effective charge defined from the Taylor-scheme ghost-gluon vertex 
in ref. ^15] (see section [2]) . We also show this low-momentum formula to describe pretty well 
some lattice ghost dressing function data [6l|13] for different volumes and /?'s. Some details 
of the computations are also provided in two appendices. 



2 The ghost propagator Dyson-Schwinger equation 



We will start by following ref. [9j and examine the Dyson-Schwinger equation for the ghost 
propagator (GPDSE) which can be written diagrammatically as 



V 



o 

k 




In addition, K-I. Kondo triggered very recently an interesting discussion about the Gribov horizon condi- 
tion and its implications on the Landau-gauge Yang- Mills infrared solutions [121113) . 



2 



that, after omitting colour indices and dividing both sides by k'^, in Landau gauge reads 



\q-kf 



k^ 



Hiiq,k) 



(1) 



where F stands for the ghost dressing function and A for the fuh gluon propagator form 
factor, 



{Al{k)A^^{-k)) = gj;,ik)5'^'^Aik') 



(2) 



with k'^g^y{q) = k'^g^v — k^ky. The standard tensor decomposition for the ghost-gluon vertex. 



n'"^{-qMq-k) 



igQf^^qu'^u'u{-q,k;q - k) 



igor""" ( quHi{q, k) + {q- k),H2{q, k) 



(3) 



is apphed, where q and k are respectively the outgoing and incoming ghost momenta and 
50 is the bare coupling constant. It should be noticed that, because of the transversality 
condition, H2 defined in eq. ([3]) does not contribute for the GPDSE in the Landau gauge. 
The integral equation ([T]) is written in terms of bare Green functions. However, this equation 
is only meaningful after the specification of some appropriate UV-cutoff A, for instance: 
F{k'^) — )• F(fc^,A). It can be cast into a renormalized form by dealing properly with UV 
divergencies, i.e. 



k^ARik^,f,^) 

Fdk^fi'') 



Z-2(^^A)5o^(A) 



Z3-^(^^A)G(A;^A) 
Z,\f^^A)F{k\A) , 



(4) 



where ^ is the renormalization momentum and Zg, Z3 and Z3 the renormalization constants 
for the coupling, the gluon and the ghost respectively. Zg is related to the ghost-gluon 
vertex renormalization constant (defined by Tr = ZiTb) through Zg = ZxiZ^j"^ Z-^~^ . Then 
Taylor's well-known non-renormalization theorem, which states that Hi{q,0) + H2{q,0) = 1 
in Landau gauge and to any perturbative order, can be invoked to conclude that Zi is finite. 
Thus, 
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Z^{^i\K)+NcZ^ gjii^i') S,^(A:^/i^;A) 



(5) 



where 



S«(A:2,^2;A) = J' 



'-<1S? j4. 



FR{q\^'')AR{{q-k)\^P) \{k-qf 



q^{q — k^ 



A;2 



Hi,R{q,k; fi^ 



(6) 



One should notice that the UV cut-off. A, is still required as an upper integration bound in 
eq. ([6]) since the integral is UV-divergent, behaving as / dq"^ / q'^ {1+llas / {27t) log (q/ f^)))^^^^^'^ ■ 
In fact, this induces a cut-off dependence in 'Er that cancels against the one of Z^ in the r.h.s. 
of eq. (0), as can be easily seen by checking that ^^^(/i^, A)Er(A;2, /z^; A) approaches some 
finite limit as A ^ cc since the ghost and gluon propagator anomalous dimensions and the 
beta function verify the relation 27 + 7 + /? = 0. This is in accordance with the fact that the 
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l.h.s. does not depend on A. Then, we will apply a MOM renormalization prescription: this 
means that all the Green functions take their tree-level value at the renormalization point, 

Fnifi^fs^) = fi^Anifx^fi^) = 1 . (7) 

The subtraction point can be taken at any non-zero scale, fi'^; however we prefer it to lie on 
the UV momentum domain (to have to deal, for instance, with renormalization constants or a 
renormalized coupling, gji{fi'^), in eq. ([5]) that could be estimated from perturbation theory). 
We will also choose to renormalize the ghost-gluon vertex at the Taylor-theorem kinematics 
(i.e., a vanishing incoming ghost momentum), thus 

Zi{fi^){Hi{q,0)+H2{q,0)) = I. (8) 

^ V ' 

1 

Now, in the following, Hi{q,k) will be approximated by a constant with respect to both 
momenta and our MOM prescription implies thus 

H,,R{k,q;fi^) = ZiHi{k,q) = Hi, (9) 

where Hi is the assumed-to-be constant bare ghost-gluon vertex. 

Although we cannot forget that the UV cut-ofF dependences in both sides of eq. ([5]) match 
only in virtue of the previously mentionned relation between the ghost and gluon propagator 
anomalous dimension and the beta function, in order not to have to deal with the UV cut-ofF, 
we can procceed as follows: we consider eq. ([5|) for two different scales, k and p, such that 
p^ — k"^ = 6'^k'^ {6 being an extra parameter that, for the sake of simplicity, will be taken to 
be small enough as to expand on it around 0) and subtract them 



Fnik^f,^) FK(p2,//2) 



9ij(Ai ) ^R{k , /i ; oo) - T.r{p , ^ ; oo) 



(10) 



Then, the subtraction renders UV-safe the integral in the r.h.s. and the limit A — )■ oo can be 
explicitely taken. 



V /■J,2 2 ^ V ^ 2 2 N rr f ^^1 i / V" , 

Sij(A: ,/i ;oo) - Sr(p ;oo) = Hi I -7—77 \ 5 — q 



(27r)4 I g2 y 
A((g-A:)2,/i2) i:^ [{q - pf , 



{q-kY {q-pf 



:ii) 



An accurate analysis of eq. (jlOp requires [5j, in addition, to cut the integration domain of 
eq. ([IT]) into two pieces by introducing some new scale q^ {qo is a momentum scale below 
which the IR behaviour is a good approximation for both ghost and gluon), 



^R{k',f,';oo)-^R{p',f,';oo) = Hi l^hnik') + Ivy{k')) (12) 

where Iir, represents the integral in eq. ([TT]) over q^ < q^ and luv over q'^ > q^. We only 
wrote explicitly the dependence on k"^ for the r.h.s. because we shall expand on 5 around 
zero with /i^ kept fixed. Then, for k'^,p^ <^ gg, we will propose the following ansat^: 

. - 4.^o(i,)) . (13, 



■^This is the massive gluon propagator where the gluon running mass [16], M(g^), appears to be approxi- 
mated by its frozen value at vanishing momentum, Af(0). 
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for a massive gluon propagator that implies of course ac = 1, as the current lattice data seems 
to point to. At this stage, we should remember that eq. (jl3p . or any other additional hy- 
pothesis about the low-momentum gluon behaviour, is needed to specify the 0(g^)-correction 
(next-to-leading) for the ghost dressing function. It should be also noted that, provided that 
the gluon propagator is to be multiplicatively renormalized, the mass scale, M, in eq. (|13p 
does not depend on renormalization scale, /i^. Now, we shall look for the ghost dressing func- 
tion, i^iR, its leading behaviour being parameterized through a general power law behaviour. 



OLF 



1 + 



(14) 



where ap > —2 to keep the integral /ir infrared convergent. After some algebra (see appendix 
, we obtain: 



^ rqo 

(15) 



where 



(g2 + A;2 + M2)2»+l (g2 + A;2)2i+1 



+ k' 



2i + l 



2i + 1 



(g2 + k^)2i+2 (^2 + ^2 + M2)2i+2 



and 



Ci 



127r24' 



r(-3/2 - i)r(l/2 - i)T{5 + 2i) 



(16) 



(17) 



Prom now on, we will focus on the decoupling case: ap = 0. Then, as shown in appendix [Bj 
the integral in eq. ()15p can be written as a series in powers of k"^, the leading term given by 



Im{k^ 



2 Aif,^)B{f,' 
647r2 



M2 



M2 



(18) 



Then, the first correction to the leading constant term for the ghost dressing function should 
be 



such that 



F{k^,fi^) F{p-^,^i^) ^(/x2) M2 
the eq. ([10]) being satisfied when: 



g2 11 



k' 



+ 



— In—-- +0[—,5^ 



M2 6 



M4' 



(19) 



(20) 



NcgRifi )Hi ^^-2 



Thus, up to corrections of the order of k^ /M^, one shall have: 

NcHiR 2 



167r 



11 



(21) 



(22) 
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where: 



R = ^F,i(0,/.^)A:K(0,^^) = ^ (23) 



It worth pointing that, provided that is renormahzed in the Taylor scheme (the incom- 
ing ghost momentum vanishing in the renormaUzation point) ii is a /i-independent (RGI) 
quantitjl^, as it is manifest from eq. (|23|) . where ax = g'|^/(47r) is the perturbative strong 
couphng defined in this Taylor scheme [17j. However, for phenomenological purposes a cou- 
pling vanishing at zero- momentum is not convenient and, instead of that, a non-perturbative 
effective charge is defined from the gluon propagator in ref. [18], within the framework of 
the pinching technique [T9] , which can be appropriatedly extended to the Taylor ghost-gluon 
coupling [15j. As a consequence of the appropriate amputation of a massive gluon propagator, 
where the gluon mass scale is the same Rl-invariant mass scale appearing in eq. (jlSp . this 
Taylor effective charge is frozen at low-momentum and gives a non- vanishing zero-momentum 
value [15] . 

aT(0) = lim (g2 + M^) ^4?^ =M^R, (24) 
ij->-o 

in terms of which the ghost-dressing-function subleading correction can be expressed: 



+ <^\W)) (25) 



It should be also noted that eqs. (|19|22|) imply to take M^/^q <C 1, as it is manifest from 
eq. (jlSp . However, any correction to that approximation will not play at the order of the 
coefficient eqs. (j23l24p . that will keep the same value disregarding that of M^/^q, but at the 
order of the gluon mass, M^, inside the logarithm (exactly like the factor 5/6 in eq. (jlSp). 



3 Comparison with ghost propagator lattice data 



In the last few years, many works have been devoted, at least partially, to the computation 
of the ghost propagator by using lattice simulations. In ref. [13j, some of those ghost propa- 
gators results were collected, mainly the ones for big lattice volumes from ref. [6], and studied 
in a different context but shown to verify the asymptotic low-momentum expansion for the 
ghost propagator where only the leading term, g^log(g^), was kept (see fig. 3 of ref. |13j). 

Now, we will consider among the results collected in ref. [13] those for the bigger lattice 
volumes and confront them to eq. (j22p or eq. (|25p . where the C'(g^)-corrections have been 
incorporated. Thus, the low-momentum behaviour of the ghost dressing function being deter- 
mined by the gluon mass, M, and the zero-momentum effective charge, aT(0), they could be 
obtained by fitting eq. (f22]l to the lattice data and applying eq. ([23|) . However, as previously 
pointed, R = aT{0)/M'^ is a RGE-quantity that can be directly obtained from bare lattice 
ghost and gluon propagators. The latter is done in ref. [13], in particular at /? = 5.7 for a 80^ 
lattice (precisely exploiting the lattice data from ref. [6]) and R ~ 10. GeV~^ is obtained. 
Thus, we will fit eq. (f25]) . where we approximate Hi = 1 and take ar(0)/M^ = 10 GeV~^, 
to the ghost propagator lattice data and obtain the curve plotted in fig. [Tj for the best-fit 
parameters given in tab. [TJ 

Thus, the ghost propagator lattice data behave pretty well as eq. ()25p asks for with a gluon 
mass, M = 0.50(2) GeV, in the right ballpark (roughly from 400 MeV to 700 MeV) defined 

''This claim is equivalent to that of ref. [13] about the cut-off independence of the bare ghost-dressing- 
function subleading term. 
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R = aT{<d)/M'^ (GeV-^) 


M (GeV) 


ar(0) 


/5 = 5.7(804) 


10(1) 


0.50(2) 


2.5(3) 



Table 1: Best- fit parameters obtained by describing the ghost dressing function lattice data with 
eq. (pS)) (see fig.[I]). The errors quoted do not account for any systematical uncertainty. 

by phenomenological tests [20] or direct lattice measurements from the gluon propagator [21] . 
It should be emphasized that, the RGI quantity R being determined by vanishing-momentum 
ghost and gluon propagators [13], the only parameter controlling the functional behaviour of 
the ghost propagator to be fitted is the gluon mass. However, it should be remembered that 
0(M^/gQ)-corrections in eq. (fT8]l will play at the order of 0{q'^) in eq. ([25]) and, although not 
modifying the low-momentum functional behaviour, the fitted gluon mass can be borrowing 
something from these corrections. Consequently, the latter prevents us to take that fitted 
gluon mass as a precise determination but as an approximative value that indeed appears to 
be in the very right ballpark. 




1.01 0.1 1 



k (GeV) 

Figure 1: Ghost dressing function from lattice data [MTB] pretty well described by the low-momentum 
formula, eq. with a i? = 10(1) GeV"^ ^^^^^ ^nd a best-fit for M = 0.50(2) GeV. 

In summary, the low-momentum ghost propagator dressing function computed from big- 
volume lattices can be very well described by the asymptotical formula eq. (|25|) with a value 
for the vanishing-momentum effective charge, aj'(O) ~ 2.5, pretty well in agreement with its 
direct lattice determinations as the gluon mass, M ~ 0.5 GeV, lies on its phenomenological 
range. 

4 Conclusions 

The ghost propagator DSE, with the only assumption of taking Hi{q, k) from the ghost- 
gluon vertex in eq. ([3]) to be constant in the infrared domain of q, can be exploited to look 
into the low-momentum behaviour of the ghost propagator. The two classes of solutions 
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named "decoupling" and "scaling" can be indentified and shown to depend on whether the 
ghost dressing function achieves a finite non-zero constant {ap = 0) at vanishing momentum 
or not (ap 7^ 0). In accordance with the fact that the lattice simulations indicates that the 
gluon propagator is finite and non-vanishing at zero momentum, we applied in this paper a 
model with a massive gluon for the infrared gluon propagator to obtain the low-momentum 
behaviour of the ghost propagator. We focussed on ap = (decoupling) and derive an 
asymptotic expression reliable up to 0{q^) for the low-momentum ghost propagator. This 
low-momentum behaviour results to be regulated by the gluon propagator mass and by a 
regularization-independent dimensionless quantity that appears to be the eff'ective charge de- 
fined from the Taylor-scheme ghost-gluon vertex at zero momentum. Finally, this asymptotic 
expression is also proven to fit pretty well the low-momentum ghost propagator data obtained 
from very big lattices simulations with a gluon mass, M ~ 500 MeV, that appears to lie on 
the right ballpark. 

Acknowledgements: JRQ and MEG acknowledge support from the research projects 
P07FQM02962 funded by "Junta de Andalucia", FPA2009-10773 and Consohder-Ingenio 
CSD2007-00042 funded by the Spanish MICINN. 

A The integral /ir 

The integral Iir is defined as: 
/iR(fe') = / 



(27r)4 I g2 y f,2 



AiR((g-fc)2,M') AiR((g-p)2,/.^ 



2^ 



{q-kY {q-pf 
where, as shown in eqs. ()13|14p . we take: 



(26) 



AIR('Z^^^) - (27) 
FiR(g2,/.2) ^ A{p''){^^Y . (28) 
Provided that ap > —2, the integral Iir shall be infrared convergent and one then obtains: 



1 1 



g2 + fc2 - 2kq cos 9 q^ + k'^ + M'^ - 2kq cos 6 
1 1 



g,2 _|_ p2 _ 2kq cos 9 q^ + + — 2kq cos 9 
1 2A(^2)5(^2) p r 



(29) 



M2(i+°f) (27r)3 

1 1 



^dq / sin^^d^ 
Jo 



q^ + k'^ — 2kq cos 9 q^ + — 2pq cos 9 

^ + „2 , ^2 , 2, o^..„.. ) ■ (30) 



g2 + fc2 + M2 - 2kq cos 9 q^ + p^ + M'^ - 2pq cos 
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We now apply that: 



sin^ e de 



+ k'^ — 2kq cos 6 + k"^ ^ \q^ + k"^ 



c. 



where we have taken into account that angular integral vanishes for odd powers of the cos. 
The coefficients Ci can be analytically obtained: 



Ci 



Jo 



de sin^^cos^^e 



127rM 



2/li 



r(-3/2 - i)r(l/2 - i)r(5 + 2i) 



(32) 



Then, one can write: 

1 2A{^l^)B{^i^' °° 



/ir(A;^) 



M2+2"F (27r)3 



j;(4A;2)*C, / q' 
i=0 -^0 



1 



2\i 



{1 + 6') 



+ 



1 + 



(33) 



(g2 + A:2)2m (^2+p2)2m (g2+p2 + ^2)2m (g2 + ^2 + ^2)2i+l 



R 



and thus expand in terms of 5, 
R = 6^ 



(g2 + A:2 + M2)2i+1 (g2 + }.2yi+l 



2i + l 



2i + l 



(g2 + (g2 + ^2 + ^2)2i 



+2 



+ 0(5^) 



(34) 



to obtain eq. (|T5|) . 



B The case ai;^ = 



When = 0, the integral /jr in eq. (jlSp can be expanded as a series on powers of fc2 
leaded by the following term: 

g(4fe2)^ Q g3+2^^^ (^^^2 + ^+^2)2^+1 - (^2 + p)2*+l 

, (2i + l)fc^ _ (2i + l)fc^ \ 4^ 

(g2 + A;2)2i+2 (g2 + ^2 + M2)2i+2 ^ + ^ 

^, ,2^(^2)5(^2) A:2 



(27r)3 M2 

/"CxD 

V 4' Ci / t 



oo 

1+j 



,=0 - v(l + i + ^)2^+i (l + t)2^+i 

2i + 1 2i + 1 

■(l + t)2^+2 " (l+i + i^)2i+2 



+7^^^^ - . . . A.2.o..,o + • • • (36) 
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Then, by integrating in eq. ()36p and expanding consistently in terms of k"^ /M'^ , one obtains: 



/iR(fe') 



,2 A{^?)B{^,^) e 



647r2 



M2 



In—,- - + O ^ 



AP 6 



2\ n 



+ 0{^^,6') , (37) 



which is the eq. (jlSp that gives the result for /jr in the regular case. 
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